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PREFACE c 


"Underway on meclear power." ‘With this simple but historically 
significant message, the U.S.S. Nautilus heralded a new age in power 
generation. Powered by the first mobile nuclear reactor, she set the 
stage for dramatic developments in military and industrial fields. A 
powerful vessel, but unwieldy and undesirably large for a submarine, 
her size was largely determined by the size of the power plant com- 
ponents. One of the more significant of these components is the pres-~ 
surizer for the reactor primary coolant system. This tank is extremely 
large due to the conservative design necessitated vy ignorance of the 
thermodynamic transient behavior in its pressure range. 

The objective of this thesis is to produce design data for steam 
pressurizing systems, by electronic analog simulation of the thermo- 
dynamic transients which occur in the pressurizer vessel. 

The writers wish to thank Professor Eugene E. Drucker, and Profes- 
sor Hugo M. Martinez of the U.S. Naval Postgraduate School for their 
assistance and enthusiastic cooperation in the prosecution of this 


project, ° 
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CHAPTER I 


INTRODUCTION 


The prototype (Mark 1) of the Submarine Thermal Reactor (STR) was 
the first full scale power reactor in the world to be completed and 
successfully operated. The first central station nuclear power plant 
in the United States will be the Pressurized Water Reactor (PWR), 
nearing completion at Shippingport, Pennsylvania. The many highly 
desirable aspects of the pressurized light water reactor cause it to be 
one of the most promising types, to date. 

A basic requirement for these light water reactors is that a very 
high pressure is maintained on the primary coolant. Pressures in the 
neighborhood of 2000 psia and higher are presently in use. High pres- 
sures permit the use of sub-cooled water at high temperatures in the 
reactor without danger of boiling. 

when any control program other than constant average primary loop 
temperature is used, a change in the volume ot the primary coolant is 
to be expected for a corresponding change in power. A volumetric tran- 
sient is consequently induced and a surge tank is required in the 
system. Since the primary coolant must be maintained at a high pressure 
the surge vessel is also used to perform this function. 

The pressurizer must be strong enough to maintain steam and water 
in equilibrium at high pressures. l1t must be large enough to absorb 
normal and accidental surges without permitting excessive pressures 
in the primary loop. Altnough thermally insulated, it must have internal 


heaters, capable of maintaining the water and steam at the saturation 


temperature. These heaters mist have additional heating capacity cap- 
able of generating steam at sufficient rate so as to prevent excessive 
pressure drop in the primary coolant during negative (out) surges. 

Specifically, in this thesis we will assume a primary loop working 
pressure ot 2000 psia, normal volume surges of 4 cubic feet and accidental 
surges of 7 cubic feet as representative values. At steady state water 
and steam are maintained at saturation temperature (636 °F) in the tank, 
but only the pressure (2000 psia) is transmitted to the cooler primary 
loop via a relatively long stand-pipe. Assuming that the installed 
heaters adequately compensate for negative surges, this type surge is 
not considered a design limitation. Therefore only those surges caused 
by an increase in primary loop volume, or positive surges, will be 
considered here. 

When designing or sizing a pressurizer, limits must be imposed on 
the various properties of the fluids in the tank. The maximum allowable 
pressure change must be determined from consideration of the strength 
of the members involved. The expected change in average density and 
therefore volume of the primary coolant must be computed. The duration 
of the surge is determined from ehanges in power, either accidental or 
deliberate. The one remaining variable, the size of the tank must now 
be determined on the basis of the aforementioned parameters. 

The pressure change caused by a positive surge may be minimized by 
introducing a portion of the surge water into the top of the tank as a 
finely divided spray. The most conservative basis for pressurizer 
design is to assume that no spray is employed, the steam is dry saturated, 


2 


and undergoes isentropic compression. Should the steam not be dry, the 
resulting final pressure would be lower than in the case of dry steam. 

The optimum and tnerefore minimum size required can only be positively 

determined after a complete understanding of the thermodynamic and heat 
transfer mechanisms involved. 

In this project, it is felt that valuable information is available, 
both in the fields of design limitation and in determining the charact- 
eristics of high pressure saturated steam. A Boeing Electronic Analog 
Computer witn associated 4-channel Sanborn Recorder (see Fig. 1) is 
utilized to simulate the pressurizer, with the above objectives in 
view. 

An important assumption which is made at the outset is that the 
driving function, specific volume (which is directly related to tank 
level) follows a cosine curve. Analysis of actual tank transients 
indicates that this aporoximation is a good one for a large portion of 
them. This type of function can be easily produced on the analog com- 
puter, while production of more exact functions is difficult. 

In the analysis of results, every attempt is made to utilize the 
principle of geometric similarity, to more nearly generalize these 


results. 
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CHAPTER J1 


DEVELOPLEMT OF THEORETICAL B{UATLIONS 


The initial step in the analog simulation of this problem consists 
of a descrintion of the system in terms of theoretical equations. 
whenever possible known thermodynamic relationsnips are used, and when 
this is not possible empirical relationships are developed. No reference 
can be found in the literature treating the following development of 
descriptive equations, and it is believed that this approach may be a 
new one. 

4S an initial simplifying asswnption, an equivalent system is devised 


which lends itself more readily to analysis than does the actual tank. 


SSSA 
N 


HEAT SINK 
2 ge |. \ 
= STEAM 
dl Peet 
Bye = eo 


Equivalent System 


Figure 2 


The actual tank is well insulated, therefore the concept of perfect 
insulation for the equivalent arrangement does not introduce significant 
error. The thermodynamic system is defined as the mass of steam in the 
tank. In the actual tank there is installed internally a spray and 
degasifier assembly. In this analysis of positive surges, without spray, 
its only effect is to function as a heat sink. The heat sink in the 
aii allem arrangement represents the tank walls, spray and degasifier 
assembly, the mass of water in the tank, and all other apparatus within 
the tank insulation. The piston represents the water surface, which 
in the actual tank performs’ work upon the steam during an in-surge. 
Thermodynamic equilibrium is assumed at all times. 


Equation #1 is the well-known work equation: 


-[r dv’ 


We 
where w - work (BTU/1b), positive when work is done on the system 
P = system pressure (psia) 
v = specific volume of steam (cubic feet/lb) 


Considering units, 
WwW = “14h P dv = -,10853f7 dv 
J 
In Heaviside operator forn, 


w- “1(.1853 P dv) ats -.1853fp (dv) dt 
p dt dt 


Equation #2 is the First Law of Thermodynamics: 
Oe= qatw 

or @= QtwWwt @& 

where e m specific internal energy (BTU/1b) 


heat flow (BTU/lb), positive when into the system 


xe) 
i 


Equation #3 is the equation of state for a non-ideal gas: 


where 4 is defined as Py (dimensionless compressibility factor) 
ciT 
T = absolute tenperature (°) 
3 
ce ee 99) (linet (sas constant) 
1 = bbb idg ft & 
18,0Lx1 44 in™ 1b_°R 


Equation #4 relates P to e and v empirically from steam table 
values (see Appendix 1): 


Pees sien? 590 v= = 30,700 v —aallen 


Equation #5 relates Z to P and v empirically, in the saturated 
and superheated region, trom steam table values(see Appendix 1): 
Z=21.715 x104P-9.47 v + 6.048 v - .5693 

Note: Equation #5 is not required mathematically, but is desirable 


for computer circuit simplicity. 


Equation #6 relates heat flow to temperature differences and is 
developed as follows: 
First, consider the sink as a system, 
dlcink = Msink Csink ST yin, 
quantity of heat energy (BTU) 


where “sink = 
Moir,” = mass of sink (1b) 
Toink = temperature of sink (25) 
Coink = Sink heat capacity (BTU/1b°R ) 
nearranging, 
sink ~ Te desink 


Integrating, 


eae = I “sink + CCE kc 
Msink°sink 
but Yetean = 7 “sink = 
and Ue) een = (To Jeter = Ty 
Then, 
Bua = Sis 2 ale 
Msink sink 


Now, considering the original system (the steam): 
dy = -K (T - Tsink) 
dt 
K represents the etfective therinal conductance of the system boundary, 
and will be considered a constant in the equivalent system for any given 
surge. The units of K are (BTU/sec°R). K is ordinarily the product 
hA, when each is determinable. 
Substituting for Dean? 


d = -~K (T + eI - T,) 


dt “Asink@sink 
Since q = me 4 » where ms = mass of steam (1b) 
dg= -K (T-T,) - K q 
dt Mig Msirnk¢sink 
The problem which is now encountered is that of evaluating the 
equation constants. Since accurate design data are not available, these 
constants must necessarily be evaluated in an arbitrary and approximate 
manner. 
For the subject tank, m.,/, can be taken as 3000 1b. This is 


approximated as follows: 


. 
‘ 
a. 
ce 
i 
a 
. 
. 


é 


mass of metal ~ 2000 lb. 


mass of water “2 1000 lb. 
Tovar 3000 lb. 


The heat capacity of water at the pressure and temperature range 
involved is approximately 2.0 BIU/1b°R . The heat capacity of stainless 
steel (tank wall material) is approximately .13 BTU/1b R . A weighted 


average heat capacity of the sink is computed as follows, 


ey = 2x 000s 17 x 2000 (Maids atu/1b F 
3000 


The mass of steam in the system will depend on the initial tank 
level, and saturation conditions, but rarely varies greatly from an 
average value of 175 lb. Properly, the mass of steam should be computed 
and Equation #6 modified for every run simulated, but for simplicity 
this average value is used for all surges. In all simulated surges, 
a, ee and m, are taken as constants, while the quantity K is 
adjusted to cause computed transient properties to coincide with experi- 
mental data. 

Equation #6 with constants evaluated becomes, 


Ose e2 x 10m = 7) = bes x 0K 


dt 


Keg 


In Heaviside operator notation, the equation takes the form, 


eo ye 10k (= T) + AAS oo toReK | 
P 


During all runs the assumption is made that the steam mass remains 
constant throughout the surge. Since the energy introduced as work 
during a surge is on the order of 3000 BIU, the maximum mass change 


would be in the neighborhood of 6 lb., if all of the energy were involved 


in a phase change. Tnrerefore, the percentage mass change can never be 
a significant value, especially since heat flows abundantly to the sink. 
Any phase change will significantly atfect K however, because of the 
heat transfer mechanism involved (heat of vaporization). 

The six theoretical equations as developed and in their most useful 


forms appear as follows: 


#1 Wos - 1 fe P (dv) dt 
i [tp ? (ey | 
if2 Oe IS i aaa 
al , 
fly Pee (2G - co V + Cp vo = Ce 
itd Zam c¢P + Cy V = Cg ye - Co 
it6 ee) * en 2] 
wie se, 
Co = Seo oe = Noi | 
hs 
c, - 3.07 x 10 Co = 5693 
E 4 a = 
oy = 6209 2a LO C10 = See wel ; 
eg = 1.715 x or oe les 
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CHAPTzR Lil 


APPLICATION CF sw U..TICNS TO THE ANALOG COMPUTER 


After the derivation of the six theoretical equations, a computer 
circuit must be designed to accomplish their simultaneous solution. In 
the circuit design .any things must be considered; such as, scaling 
factors, follow-up tiie of multiplier co.ponents, amplifier drift, volt- 
age magnitudes, etc. A certain amount of latitude is allowed in circuit 
component selection, but detinite limits are present due to stability 
considerations. 


The six basic equations in co:aputer notation appear as below: 
: P 


ral eee es Teer (eae, de.) 

p dt ; 
#2 S & soees a eos ses e 

ae es —\ 
#3 es 3-357U 8 | 
iG 

fh fee 255 = = 5.167 ye 6+) «6453825 (.02 v) = 5200gee 
#5 WOU 3 ones 486 - 2.3675 (02 v-) = 28.065 
ib G s-L[lmekT - 176k T + .000u5 K 4 | 

Pp 


The bar over a yuantity represents the quantity as a voltage. The 
detailed conversion of the theoretical equations to the above "machine" 
equations couprises Appendix il. 

Fig. 3 shows the resulting machine circuit in olock diagram form. 


For the detailed circuit, see Appendix li. 3riefly, the circuit is 


desicned to function as follows: 
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In the driving function circuit a cosine function is added elect- 
rically to a constant voltage, the resultant voltage transient very 
closely approximates the specific volume throughout a surge. The 
amplitude of the cosine function represents one-half the change inv , 
while the period of the function is twice the time of the surge (4+). 
Within this same circuit the time derivative of specific volume, dv 
is also generated. = 

A function multiplier is used to obtain the quantity Y, while a 
servo-multiplier produces the products P dv and P v. An integrating 
circuit integrates P dv as a function cee to produce specific work, 
w. Another ieee circuit produces the time integrated specific 


heat flow, q , using temperature inputs. 


A summing amplifier sums a, wy, and q to produce e. Another sums 
ve 


e, v, and v~ to produce P. A third sums 2, vy, and to produce ae 


A dividing circuit performs the division to produce T. 


Py 
Z 
The circuit required to solve the set of equations consists of 
eleven amplifiers, four sign-changing amplifiers, two function multipliers, 
one duo-channel servo-multiplier, and twenty-one 50,000 ohm potentio- 
meters, The net result is a rather complex circuit, but one which does 
not seem to be subject to further simplification without sacrificing 
necessary accuracy. Since instability in this type of computer is a direct 
function of the number of components used, circuit simplicity has been 
a constant objective. 
This phase, the designing and setting up of the circuit, requires 


a great deal of computer experience and proficiency. For the uninitiated, 


many trial and error situations are encountered, and with a circuit of 
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this complexity the process can be very time consuming. 


The following chapter will deal with the problem of producing 
quantitive results, assuming the circuit has been designed, balanced 


and tested. 


CHAPTER 1V 


ANALOG COkPUTER SOLUTION OF THEORETICAL KLQUATIONS 


After the analog computer circuit is designed, assembled, and rough 
qualitative results obtained, the next phase consists of producing 
accurate quantitative solutions. The first step in accomplishing this, 
consists of static adjustment. With various initial values of ey» and 
v set, the e, P, Z and T computing circuits are adjusted to give steam 
table values for these quantities over the entire anticipated range. 
This step ofthe co:puter set-up is very exacting and time consuning 
but is a necessary step prior to adjusting for dynamic accuracy. 

Prior to the dynamic adjustment, a hand solution of the system of 
theoretical equations is obtained for use as an adjustment reference. 

A set of initial conditions, and a driving function are arbitrarily 
selected, and using the Heun Method of numerical integration, a point- 
by-point solution of the equations is calculated (see Appendix II1). 
The various thermodynamic quantities as calculated are plotted, the 
resulting curves serving as the final reference for dynamic adjustment 
(see Fig. 4). 

Next, the computer is made to duplicate the hand solution. This 
is done by adjusting the w and q integrating circuits until, the w and q 
curves, as viewed on a properly calibrated Sanborn Kecorder, are of the 
proper shape and magnitude (see Fig. 5). These adjustments in general 
are minor. 

The value of careful static adjustment is now apparent since no 


further adjustment of the e, P, Z and T circuits is found necessary. 


1 Special case of Runge-Kutta Method; see Kopal, Numerical Analysis, 
pp. 202~205, Wiley (1955). Lh 
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The curves of these quantities are now round to match tite correspondin: 
hand solution curves very well (see Figs. 4 and 5). 

Once the hand solution is duplicated, a reference tor measuring K 
is now available. In the hand solution, the arbitrarity selected value 


of K is .9254. The corresponding setting of potentiometer a,, is now 


14 
readily observed (usually in the neighbori.ood of .100, varying slightly 
from day to day). Since K and the ary setting are proportional, K can 


be determined at any time by the simple relation, 


ho Kh. x (ayy) 
Oo lem 


It is now possible to vary K over a larse range of values by merely 
adjusting a single potentiometer, thus varying the intecrated , curve, 
and all other curves of thermodynamic properties. By sinulating the 
driving function for a given surge, K can be varied until the proper 
pressure curve is obtained (i.e. a oressure curve which matches the 
pressure curve resulting from an actual tank transient). The neat 
chapter discusses the investigation of actual surges by si.uwlator, to 
determine Hallie of K for tnese surges. Correlation of K with surge 


types is also discussed. 
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CHAPTin V 


LAPSRIMENTAL RESULTS AND CCNCLUSGICNS 


Appendix IV is a compilation of data obtained from the actual surge 
tank. wighteen positive surges are represented, spray being utilized 
in only one. It is evident tnat possession of many more runs would be 
desirable, but the present number must surfice since ore could not be 
obtained. fortunately these eighteen surges represent a fair cross-— 
section of positive surges encountered. 

The procedure for simulating a given surge is outlined in Chapter 
lV and briefly is as rollows: 

A proper driving function (specific volume) is set up on the comput- 
er, providing for the generation of e cosine curve of proper amplitude 
and period. Initial values of v, e, T, and a trial value of K are set, 
and a surge is generated. K is adjusted for several trials until the 
generated property curves match the corresponding actual tank curves 
(see Fig. 6 for Sanborn recording of this procedure). 

Each of the eighteen surges is simulated in this manner, and the 
proper value of K-~determined (see Table 1). many attemnts at correlat- 
ing K with various quantities (av, av , av , etc.) were fruitless, but 
a fair correlation with surge Patio a) is possible (see rig. 7). 
Mis lixely, in fact, that k is a function of several or all of the 
variables considered, but the effect of each seems small in comparison 
to the effect of at. By Siuple curve fittinc, a correlating equation 


mo cCerived, 


h = B75 + 2 
tee Ly 
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Nothing is claimed for this equation except that it describes a K which 
aids in predicting pressure rises on the safe side (i.e. in most cases 
the oredicted pressure rise will be somewhat larger than the actual 
pressure rise). Column 4 of Table 1 lists the K for each surge as 
predicted by formula. 

in general, it can be observed that K does not vary greatly for 
any surge of over 30 seconds duration, and that it becomes nearly cons- 
tant for surges longer than 90 seconds. it is also true that pressure 
rise is not sensitive to moderate changes in K. The last colunn of 
Table 1 lists errors encountered when a constant value of 3.75, the 
limiting value of K, is used in the computer. 

Tre mechanisn of heat transfer for the actual tank can only be 
surmised. However, it is logical to expect greater turbulence for short- 
(©. iasver surzes, resulting in somewhatm@iarcer values for kh. This 
could be due to a larger leat transfer suriace between the stean and 
water in the tank, The effectiveness of spray in increasing K is readily 
observed in run 1506, Table 1. Whereas a Kk of approximately 5.0 could 
have been expected without spray, the actual « observed is 13.8 

From the general characteristics of the K vs. Surge Duration curve 
(Fig. 7), it appears that K approaches a limiting value of about 3.75 
for surges of long duration. For geometrically similar tanks utilizing 
the same pressure range, the same limiting value of K could reasonably 
be expected. 

Realizing that a correlation exists between 4K and the driving 
function, it is now possible to predict pressure rises for all conceiv- 


able driving functions. This is done in the following manner: 
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# \srror in predicting pressure rise (in psia) when usi:g a constant 
value of K,( the liniting value cof 3.75), rather tnau the exper- 


dental or formula values. 
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A driving function is eitablisned, the pro.er value of K is set, 
aud the resulting pressure rive recorded. in an attempt to use tie 
principle of ceouetric si.ilarity, dimensionless paraweters are used 
mien possible. sAccocdingly, a ta“ily of curves is obtained relating 


my ,AP , andat (see Fig. &), 


v E 
ae 0 : 
Since the actual tank surses vary fro. 27.5 secords to 120 seconds, 


the predictions are only valid within this ranse. However it appears 
that safe extrapclations inay be made. An attempt to justify liaitine 
values of tne prediction curves is bxsed on the following reusoning: 

initially the syste. co.sists of dry saturated steam. «as a surge 
progresses work is done on the system, and heat will flow out provided 
a temperature difference exists. lxanination of Fig. 9 shows that the 
isentropic compression process prescribes the practical upper limit of 
Over-pressure. it is believed that a process following the saturation 
line prescribes the lower limit of over-pressure for this tank. Since 
Meessirge introduces colder water into the tank, a potential sink exists 
which could conceivably cause the process to enter the wet stean region. 
However, a careful study of actual surzes indicates that the steau 
ren2ins superheated throughout an in-surse. Very long surges ap ,ear to 
parallel closely the saturated line with slisnt superheat. i1n addition, 
all surzes apnear to ve initially isentropic in nature. 

Based on the aoove reasonirs; and actual test results, it is believed 
that the oredicticn curves approacn finite values as@t approaches in- 
finity. Tse limiting value of ai’ (for any given Ov } is that which 
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For very short surges, tne liBitime conditiwy. appe@ars to ve ome 


of isentrovic couwpression, thus establishing: a Suximun villue of A 
a 


for any sivenfiv. (n the actual installation, surses of less than 10 


< 


seconds duration «re unlivelLy aque to oripary loo» circulation tise 
considerations. Thus the extrapolated curves for values of at from 10 
seconds to 30 secords would see: reasonable. Using the prediction 
curves (Fig. $), Table 2 shows a cowparison of predicteu and actual 
pressure rises. 

It is believed that basic objectives of this project have been 
realized; however many avenues for further study appear to be open. 
Although a procedure for rational analysis nas been indicated, it appears 
that a study of the effect of spray, and of negative surses,by these 
methods would be of value. 

In the field of heat transfer, a creat deal remains to be done. 
The combining of all heat transfer mechanisms into a sinzle coefficient 
is perhaps an over-simplification, althougn effective in this analysis. 
lt is known that several thernaal paths exist, and a more conplete: 
analysis could perhaps isolate the effects of each. Heat can flow from 
steam to metal, from stean to water (directly or by coidensation), and 
from metal to water. A comprehensive sinulator study could provide 


Knowledge in this relatively unexplored field. 


27 


Page. 2 


— 
RUN WO. : PLUMES Ut Oe | UG: PREDICTIUN shRut 


— te S eiaeied 
a a ne ne a a 


2302 


2045 
P28 
+1556 
2229 
1506 


2023 
2010 
2248 
2210 
2247 
2113 
hol 
2255 


PRINCIPLES UF sNGINELRUNG THESACDYNALACS 
2nd Ldition, #iley, 1954 


neenan, J.'l., and PHORMCDYRALIC PROPER oS OF GIEAL, 
eyes, F.c, Wiley, 1936 


dakeb, li., and ELEMENTS UF HEAT TdAncFun AND INSULATIONS 
Hawkins, G.a. 2nd Bdition, jakey, 1952 


Beecidans, i. HAT DiNteatoolcn, Jud wedition, 
neGraw-Fill, 1954 


Bermeoler, n.C.H, BASiC THsORY CF Tits cle T.u0iiuC ANALUG 
CU.rUTIgepenner oclertarie Co,, 1955 


APPENDLA J 
DERIVATAGN CF EMPIRICAL RELATIUNG FOR PiteSOURE 
AND COwPheoSIBILITY FACTOR 
Since e can be computed from the First Law, e = gq + Wt 8 
and v is available as a driving function, it is convenient to express 
P as a function of e and v.. This can be done by constructing an 
empirical relationsrip from steam table values, which will be accurate 
within certain limits. The pressure ranze which is considered sigznifi- 
cant for this problem is 1700 - 2700 psia, and no accuracy is claimed 
beyond these limits. 

Since e cannot be taken directly from the steam tables in the 
superheat rezion, a plot of enthalpy vs. pressure is made initially, 
showine lines of constant temperature and constant superheat (see Fig. 
10). 

hext, at various te.peratures, various values of h,pe , andv 
are listed. For each of these values e is computed from the relation 


& 


= h - 144 Pv (see Table 3). From this tasle, a new chart is 
constructed <aieate e,P,v andT (see Fig. 11). Using this chart 
as a basis, initially an expression for e =f (P,v) is obtained in 
the following manner: 

The average slope of constant v lines is .1176. For any one line, 
meee .lL76 P + Cl, The constant Ci), is deterinined as a function of 
v (see Table 4 and Fig. 12), = 562% 20536. The resulting 
expression for internal energy, ee i76p + 1562 v + 538 , when 


tested at the extreme limits is found to sive considerable error. 


Although undesirable from the standpoint of computer circuitry, au 
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FIGURE 


TABLE 3 


P (psia) v (£t7/1b) e (3IU/1b) 


—_ 


62553 al 1105.6 

sau vee? L800 £2296 Tbo.6  ' 1Ooimma 

1145.8 L400 «2040 lees | luvaae 

(sat.) ye 1919.3 .1992 10.3 | 10703 
6,3 USS 1800 mO HO BO.3 | LlOSme 
(Oe 1900 Eolas on Loyoee 

ESE 2000 1936 aie | 1073s. 

(set.) 1130.6 2059.7 1798 68.7 1061.8 
650 1167. | 2000 2058 (6.2 1090.8 
1146.3 2100 .1846 False 1074.5 

TALS 2200 1633 6.4 1054.6 

(sat.) 11.5 P2082 ~1616 66.2 105244 
660 SLSR Y ICC 1962 1on3 LO9L 
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670 150.0 | Pektere | 702 72.5 LOT dee 
1123.2 21,00 £1526 67.8 1060.4 

1099.8 2500 lee 62.2 1037.6 
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650 ieee 2500 ue 68.7 106256 
IMEC rel 260 L539 SOS, LOlSee 
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(sat, ) 1067.2 2708.1 El 56.u LCi 
690 1156.6 250U 0159 (Oe: 1082.8 
‘eye 2600 Sly hy: 65.8 Worse) 

MLS3 200 .1299 61,68 LUAI.D 

700 1160.6 2600 1549 74.6 1086.U 
12.5 2700 nS 70,7 107 
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additional modification is necessary, making e = p(P,v, v*). 


To accomplish this, the expression P - &6.5le = C15 y? + C7 
+ C4 is set up. Substituting actual steam table values at three 
points ( P = 1900, 2100, 2300) three simultaneous equations in Cis ; 
Cy~e » and C14 result. Solving for these constants and substituting, 


the final equation for P is : 


Be Gs5i.e) + 525500 See = 161 


The calculation of an empirical relationship, 24 = f (P, v) follows 
Similar lines. First Z vs. P is plotted from the steam tables, 
showing Lines of constant v (see Fig. 13). The average slope of the 
constant v lines is 1.715 x 10-", the equation for any one line being 
meee (Lf Xx lens i ear Clg Cig is determined as a function of 
= 2.943 v - .318 (see Table 5 and Fig. 14). The expression 


1S) 
4 Pee 94 y, 


mor 4 as a function of P and v is, Hx dh TALGY ee 105 
= septs; 

Once again a v term is found necessary fcr desired accuracy, and 
using the same procedure as before the resulting expression for Z is: 


Ze = 1. fd x No Pos OOhy a Sonos Vv we «9693 
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APPENDIX 11 


DurkLVATION CF MACHING b.UnTLONS 


Sealing factors are assigned as follows: 


=) Pp - = 
x5 = saan - “=. - 60 
Pinax 
ad C 
pA Sich a ae ee — er 1 
oe 50 
max 
a7 = Trax = 1500 = 30 
r.  » 
Max 


A 50 
max 
PS ceeeeemer = 1 = .02 
A 50 
ek 


Wmax 
max 
x = if = ab = ile ' 
v t AL 


ie | -2851 ‘e P av | at] 1O o, 0, ¢ 
ne ae wt 
ee: a [ petut 1 P av) at | 
p dt 
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eau fe Cava cg ¥ - 
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p 2 


culated potentiometer settings, and components selected: 


- OD. F. dependent ag = A255 25 = .4383 
= iN " ag = 6864 a16 = - 6864 
a ‘ Big a aoe Qj = +5000 
= .llil aj, = .6333 @jg = +0089 
Bees 333 aya = 302K Stig =e 

= -3357 ay3 = .51K5 2 = "2368 
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HAND SOLUTIUN UF THsORETICAL L YUATICNS 


The first step in obtaining a hand solution consists of combining 
the six theoretical equations to form two dependent differential 
equations: 

Substituting from equation #2 into equation //4, and then into 


equation #1 one obtains, 


du [ AL 
SFG gh eo Cay os G3 os os [dv 
dt. ee) * Seu. 2 | ae 
Since Va, te Se wt 
dy 
en de = ~Usw Sin wt 


and coz: ; 
Gi Uowsnat CC, [qgrw-¢ fy + Stu's e,- Se | 
un as qe 
= z 
Substituting from equations /4 and #5 into equaticn #3, the 
resulting expression for T becomes, 


ee ce ae 
+s pelo Fo ev+ Gy + @, oo! 


ae ae ~ oe = =) Zo = & | 
[3 1, - ev t(a va a - 4.) 
it is now convenient to ae values to the Fin. tunction, 


initial conditions, and to evaluate the constants: 


uy = .OL9 ft/lb cy = 8.51 c, = 6.048 
Ww = .0418 rad/sec Cy Eemmro.U7 xX iy4 Cg an Jang 
Vnean = -172 it?/db C, = 5.299 X Lo# Cy = 9693 
Cg = 1066 BIU/1b co = 316L Cy> = -1853 
ey = «5957 Cae ie 107" 


Evaluating, 


d ToT 


“4 ome a as 
oat Seco 5 HO int gig t 9 +W- S3G6iI0V + GI6V ve 694] 


es 7800 v E +W- 3610V +G6iB0v +694] 
f qGtW +S4ov ~ 310 v~ + 304] 


Equation ,/6 is, 
| fe 
pe Co KT + CUR TC ia 


Arbitrarily assigning a value to h , evaluating constants, and 


substituting the already derivea expression for T into equation #6, 


K = eo 
' a) 
Sim = Sa 72 « LO 
=e 
Cul = [so leD pa®) 
Dae eg.) Rk 


== i, if -~ 3 OV +6!80V' +694 = {- 
an ove [q +w - 3610 v releov' +694] 4.12410 g pa 
L4, rW + S40 V - 3iov™ + 304 | 


The two resulting aependent first order differential equations 


Ipecone, 


d = K, [4 +U +k, ] ~ 4.12410 'G 45.770 
Lg + UW t+ Ka] 


i ete # irivin : 0 
where K, Kz 4, hare tabulated values of VY, the driving function 


(see Table 6). 


Ad 


TABLE 6 


Ky 

| -9.9U | 

ULE “d -~9.88 

‘ : ~¥.82 
Le 9.1 | 
, : -Y.56 ! 

| 0 ; | ea 

30 UU59 : Hes | ~9 24 

SB -U02U | L7hu | -Y U2 

: 40 ~.0020 | 1700 | =2.02 
45 ~.U059 PLCGL | -2.62 | 
5U -.U095 aaledes -3.37 : 
55 eecioy =), 1593 -2.26 | 
60 ol, | 1567 : -£.13 | 
65 -.UL74 ees) : -5.U2 
70 - .UL26 Boy | ~7.99 | 
75 ~ ULYO : 15 3U | (ae | 
80 eei6 1534 | ~1.45 | 
25 : afoul 1547 | 26702 | 
90 | -.0154 e Wes7) -3,13 : 

95 Sol27 «| «4.1593 ~8,.26 
150 | - .0095 | L615 ~8.37 
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The next step consists of nuuerically intesrating the above 


equations to obtain point-by-Loint values of 4% > 4 » UW {aed do , 
ch dt 
Tne Heun lethod of numerical intezration is used, the general fori 


being: g Se oe ee 
qy f Cu, giv) 
W = g (wigiv) 
(@ 


n) (n-1) 


BF hint BUG t §: 
fy 
where OE = Bent + h (alae. 


u 


(4) ~ Cn-ty 
7 h is 
EO ey Pe ( Wei oe ) 
where fry ‘ 
CU, a We 1 h (w..) 


h is the interval between »oints and is selected as 5 seconds. 
Table 7 lists the resulting values for arih, Ww, and W 
- ax 


obtained over a tine of 100 seconds (21 points). 


Next, the resulting values of 4 and W are suvstituted back into 
the original theoretical equations, resulting in point-by-p point values 
for the various properties tnrouzghout the run. These values are tab- 
ulated (see Table &), and are plotted as a function of tine in Fig. 4, 


Chanter iV. 
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TABLE 7 


HAND SULUTLUN (kh =.9265) 


| 

ia | 2.28368 16.0765 0779 
15 75 Ei O5 71 oe] Leen 0 

16 | 80 -15.6849 wagon 16.0843 -. OFie 
1 | 85 SiGe y ih -.1536 15.5423 -. lies 
18 90 SAG allen -. 1087 Lh 6944, - 1972 
19 95 -17 5063 70246 13.6080 -.237k 
20 | 100 -17.4614 0426 e375 -.2659 


PROPERTY VALUES FRUL! HAND SOLUTION 


TABLE 8 
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APPENDIX 1V 


ACTUAL TANK TRANGLiNT DATA 


Run #2247 Mio "= 9169.06 1b. 


Time Pressure 


0 2002 31.7 11875 

5 2012 SVL | 1867 
10 2042 See 1843 
1) 2094 30.4 ei 
20 2168 29.1 1720 
25 2233 Bes 162 
30 2280 26.1 1542 
35 2305 oo 1500 
40 2316 25.1 1485 
Run #2240 m, = 167.46 1b. 
0 2002 32.8 lel 1875 
5 2007 | 33.0 31.3 1870 
10 2026 3375 31.0 11850 
15 20n7 | Bhat 30.6 1827 
20 2071 34.9 30.2 1803 
25 2091 35.6 29.7 77g 
30 aon | 36.3 24.5 .1760 
35 2112 Bors 29.0 1731 
40 2112 37.0 28.9 £1726 | 


20 


Run #2229 Mo = die el Seb. 


Sie Pressure Level Volume Spec.Vol. 


0 1992 31.5 61885 
10 2007 31.3 -\ a 1873 
20 | 2033 30.8 1843 
30 | 2070 36 1 | 1€00 
10 2107 29.0 | .1735 
50 2170 28.0 | 1675 
60 2220 BO 1615 
70 2264 41.8 26.1 1560 
80 2295 42.7 25m) | 153K 
Run #2128 Me W69705 1b. 

0 | 2000 2203 | Se | 1875 
10 | 2022 32.7 ae | 1861 
20 | 2037 33.4 31.0 | 1832 
30 2072 Bye) | 30.4 1796 
1,0 | ein) 36.0 e955 1743 
50 2162 37.8 28.5 1684 
60 220i Bo. Balk 1620 
70 22hh, 41.0 26.6 ae 
80 2286 42.0 26.0 1537 
90 2299 42.6 258 1520 

100 2299 42.7 25.6 | 1514, 


ok 


Run #2113 Hes When. «Lb. 


Time | datas Level Volume | Spec.Vol. 
@) | 2000 1878 
Zo ! 2002 870 
> | eolg 1854 
i) 205k, - 1830 

10 2090 «1795 

12.5 2124 L765 

is 2S of oo 

lee) 21f2 . 1703 

20 2182 . L680 

22.5 2186 .1662 

25 ZS -L65o 

215 2187 «1645 

30 | 2187 | - L640 

Run #1101 m, = 169.5 1b. 

0 1990 eine | 32.1 wee 
2 195 B18 32.0 1888 
2 2009 Sas) Silat .1870 
lee? 2039 eae Dine ~1854 

10 2066 3300 31.0 1630 

ae.) 2085 34.4 30.4 ie 

1 2097 Bio!) 90.1 L775 

17.5 2102 35.3 29.9 1765 

20 2102 B26) 2928 -1 760 


De 


mum 72255 Me = yo 


Time = Pressure a Level 


iL Voluwe pec Vol. 
= 
0 1995 One $3.0 1885 
2.5 7 2002 ay 30.1 33.0 1885 
5 | eral 30.3 32.9 1880 
7.5 2032 | 3.6 | oom | 1870 
10 | 2043 otha! eee .1850 
12.5 | 2055 | Bley | 32.0 1828 
15 : 209, | 32.5 ! Bie : 1805 
17.5 2121 | Bee | uae 1775 
20 | 3h 34.0 | 30.7 1755 
B25 | 21u3 BG cre .1730 
25 | 2146 35.0 | 30.0 | oe 
en. 5 | ally | ae | 29.9 | .1710 
30 | 247 | 35.5 | 29.8 1703 
Run #1506 fy = pees 
0 1952 30.8 | aeae | AGL 
5 1558 3u.9 32.5 | .1940 
10 1988 Sry 32.0 | .1908 
15 2032 Gee sie 1867 
20 2082 Pol 30.4 1412 
25 2130 35.7 29.7 L770 
30 2174 7a 28.8 | 1718 
35 ale, 38.5 28.0 | 1670 
40 | 2238 Nl |) eres: | 1629 
15 f 22h 40.0 | Pao | 1616 
53 


Run #1556 


“a 


iia = Ste =} 


s 
20 
eo 
30 


on 


Run #2010 


| 
> 
20 
25 


Mg. = aged eu. 


Pressure Level Yolune | Spec.Vol. 


1893 Bile 3208 | 2030 
1894 Sine | 523 | 2022 
1923 32.0 | 31.9 .2000 
1960 33.0 pie | 1962 
1998 Bie | 30.5 1912 
2034 3504 24.8 1868 
2063 36.5 29 3 | 1837 
2086 Cae. 28.7 1798 
2101 38.0 28.3 | 177k 
2114 Corr | 28.1 17e 
D122 39.1 27.5 | Weg 
2122 | 39.4 Bile | 1718 
M, = 64.67 1b: 
19390 | 31]. 1975 
1933 31.0 : 32.5 1975 
19h2 Slat B25) 19 
1962 31.5 Bee .1951 
1988 Be.3 31.7 1927 
2021 se | eu2 .1898 
2057 3h oly B05 61855 
2094 35.4 29.8 ete 
2124 36.4 29.3 .1782 
2146 27.0 28.9 Let 
2159 37.6 28.6 1740 
Pile? By ay 28.4 en 


ok 


Ee —E EEE 


Run #2023 


Run #1630 m= 198.3 Ib. 

0 | 1938 | 20.0 | 39.0 | £1967 
10 1952 | 20.6 | 38 LIE 
20 | 1982 | aes com .1920 
30 : 2022 | 22.8 eyes .1880 
1,0 | 207k | 24.8 sroral . 1830 
50 age. 26.7 35.0 ee 
60 218h 27.8 sinh Nie 
70 222k 25.8 Spee .1670 
60 2252 ah eee £1623 
90 2262 B2n6 31.5 61589 


wun jf1911 emp eye, ela. 


Level 


pee 


es 
co 
NC re No) costo CT 


OF. 


56 


5 ele Pee 
Run #1656 


Tine 


awe 


Me 
ie 


a 


Run #2248 Mg vw 156.3 1b. 


Pressure Volune opee . Vol 


= 


Bie) 12007 
aay eecly, 
oa -2017 
«2008 
. 2000 
1982 
1958 
odie 


Run #2302 Mm 62.5 Lb: 


Tine Pressure Level 
ne a ea pee gree toe 


0 | 18&8 


1891 
1902 
1918 
1935 
i? 
1) 
1996 
2017 
2038 
205k, 
2062 


2062 


Run #2045 Mee) 16720 Sb. 


Time | Pressure Level | Volune | Spec.Vol. 
SSS Se a ey SS 


Za | 

D 1990 Dee BIL L894 
10 2005 o2a0 31.4 1880 
a) es, 33.5 See : 1856 
20 2054 3h.d SIRI os) 
2S) 2091 3549 29.8 78> 
30 2130 SISAS Coy sk 1746 
Bo 2166 Bled 28.6 TEGO, 
4,0 Ae atl Sie) 25.0 1674 
A5 Cee 39.9 27.4 1640 
510) 2252 OES Dee . L610 
DS 221) ill 8) 26.5 154 
60 | 2290 Hl.g 2G.1 1562 
65 : 229) eee Cone Laks 
70 2295 Zie2 ub eat Be 


60 


Run #2039 a ee eo = eee 


ie) 

Time | Pressure | Level 
0 1981 | BenT 
25 1981 32.0 
5 1986 33.0 
5 1995 2302 

10 | 2008 | 33.6 

fe .5 | 2025 | halk 

i | ene | Brey 

oy. 5 | BLISS; Boel 

20 | 209 | 36,1 

ed | 2137 | 37.8 

25 Pay | 33.2 

By 5 2 | Bee 

36 | 2252 oes 

32.5 | 2286 | 41.2 

oe) | 2312 | 42.4 

Bi} | yoe2 : hoa 


_ 


£8980 
Bosley : 
Simulation of steam 
: pressurizing tank trausients by 
analog computer. 


